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Abstract. 

We investigate the effect of general time transformations on the phase 
synchronization phenomenon and the mutual information rate between pairs of nodes 
in dynamical networks. We demonstrate two important results concerning the 
invariance of both phase synchronization and the mutual information rate. Under time 
transformations phase synchronization can neither be introduced nor destroyed and the 
mutual information rate cannot be raised from zero. On the other hand, for proper 
time transformations the timing between the cycles of the coupled oscillators can be 
largely improved. Finally, we discuss the relevance of our findings for communication 
in dynamical networks. 
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1. Introduction 

Time and complex dynamics play a major role in biological, social, economical, and 
physical systems. Cycles of different periods often govern their dynamical behavior 
and determine their intrinsic activity. A variety of processes require a precise timing 
between the oscillators cycles for a proper functioning, as for example, the respiratory 
and cardiac systems [1] , spike discharges and information transmission [2l [3] in neuron 
networks, ecology [I], fireflies blinking after dark, and peacemaker cells of the human 
heart [5]. Synchronization is an efficient mechanism to generate such a timing 
[HEiEllUElEilTlIElIHl EH]- Among several types of synchronization recently 
found in complex systems [8], chaotic phase synchronization (PS) displays special 
importance because of its weak constraints on the dynamics and coupling strength. 
It has been reported that PS mediates the process of information transmission and 
collective behavior in neural and active networks [2], [Til H2L as we ^ as communication 
processes in the human brain [2[ [13[ [H] . 

In real systems PS is the most common type of synchronization [lj [2j [U [5j [8j [9j 
[TOl [TT1 [T3l [T4] . The main reason for PS to be so common relays on the fact that real 
oscillators are not identical, but have some parameter mismatch. When real coupled 
oscillators undergo a transition to PS the timing is not precise. In many situations one 
wishes to improve the timing condition, but in fact, one cannot systematically control 
the oscillator parameters to drive them to a higher level of PS. The question is then 
how to improve the timing without changing the oscillator parameters. The natural 
candidate is a time transformation. Could one enhance a better timing by changing the 
time? Or even better, could one introduce PS by time transformations? 

Coupled dynamical systems under time transformations are important in physics 
without an absolute time as well as in situations where the time cannot be directly 
obtained, as in the study of sedimental cores in the field of geophysics. In the latter 
case, the time at which the sedimentation took place is usually unknown. Only a proxy 
for the time can be derived from the measurements, which does not yield the "real" 
time but only a monotonous transformation of it [15]. In the study of synchronization 
phenomenon in such a system the natural question is whether not having access to the 
real time could effect synchronization. 

Such time transformations (typically nonlinear) have attracted a great deal of 
attention (see [TH [T7j and references therein). They cause no change in the topology 
of the dynamics, but the duration of the cycles can be drastically modified. An 
important problem is to analyze whether the dynamical properties are invariant under 
time transformations [16] . Recent results have shown that dynamical systems under time 
transformation can present nontrivial and counterintuitive properties. For example, a 
nonmixing dynamics can be converted to a mixing one p2] . 

In this work, we show that time transformations, satisfying simple conditions of 
integrability, can neither introduce nor destroy the phenomenon of PS. We also explore 
the natural connection between synchronization and information exchange in coupled 
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oscillators. We uncover the transformation law for the mutual information rate (MIR), 
the rate with which information about a node can be retrieved in another node. If the 
MIR is zero in one time frame it will remain zero for any other. On the other hand, if the 
MIR is nonzero it can be drastically modified by a time transformation . Surprisingly, 
if there is no synchronization (to any extent) between the nodes forming a network, 
time transformations containing information about a particular oscillator (node) of the 
network cannot be used to carry this information to another oscillator. 



2. Two Oscillators Case: Enhancing a Precise Timing 

We first illustrate our approach for the paradigmatic example of two coupled Rossler 
oscillators: 

±1,2= - a 1>2 yi,2 - zi,2 + e(x 2 ,i - x lj2 ), (1) 

Vl,2 = "1,2^1,2 + 0.151/1,2, (2) 

ii, 2 = 0.2 + z lt2 (x lj2 - 10) (3) 

with cti = 1, and a 2 = at\ + Aa 2 . We shall denote = (xj,yj,Zj), where j = 1,2, 
and x = (xi, X2, xi, X2). Since for these oscillators the trajectory revolves around one 
specific point [Fig. ([I]) (a)], we can simply define a phase by tan0,- = Vj/xj, which yields 

m 

<A? ( x > t)= (yjXj - Xjyj)j{x) + y])dt. (4) 







Furthermore, let us denote the time at which the oscillator x^ completes its ith cycle by 
ty That is, the times at which the phase is increased by 27r (see Appendix A for more 
details). We can show that there is PS if, and only if, we have 

\t\-i\<K : (5) 

where k is the minimum finite number that bounds the inequality. For more details 



concerning this equivalence see Appendix A The value of k shows how well paced both 



oscillators are. The smaller the value of k the better the timing between x x and x 2 . 

For e = 0.0015 and Aa 2 = 0.001, the two oscillators are in PS, which means that 
the phase difference A(f> = 4>\{t) — (f) 2 (t) is bounded for all times. Consequently, Eq. 
(J5j) holds [Fig. []Jb,c)]. In the PS regime the oscillators have the same mean frequency, 
namely 

(ii)t = (02>t ~ 1-035, 
where (-)t is the time average with respect to t. The average period is given by 

(Tj) t = 2n/(<j> j ) t « 6.067. 

We have that max\t\ — t 2 \ corresponds approximately to (3}) t /4 [Fig. He)], 
which can be rather problematic for a reliable communication system based on chaos 
synchronization, since the two oscillators do not reach the Poincare section with a precise 
timing. See Ref. [H] for a detailed discussion. 
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Figure 1. The effect of time transformation in two synchronous Rossler oscillators. 
In (a) we show the attractor projection onto the (x, y)-plane, and in (b) the time series 
x\ and X2 vs. t. In (c) we depict the time difference \t\ — t\\ for ith period between 
both oscillators. One can see that even though the quantity \t\ —t\\ is bounded, it has 
large fluctuations. In (d-f) we proceed a time transformation t — > C, given by Eq. (|7|). 
In (d) we show that attractor projection in the subspace (x,y), and the time series x\ 
and 22 vs. £ in (e), while the amount \C,{ — Q\ is shown in (f). One can see that after 
the time transformation the timing condition is drastically improved. 



The main question is then what could happen when we transform the time. Could 
we produce an effective improvement of the timing? If we linearly scale the time 
introducing ( = at, then the average period transforms as 

irp \ \ T j)t 



f l — t l 



a 

while the timing 

ici-ai 

a 

Thus, there is no effective improvement of the timing with respect to the average period, 
since 

iCi-di _ 1*1-41 

(Tj)c (T ) t ■ 

The situation can be altered to improve the timing condition in Eq. (jSJ) by using a 
nonlinear time transformation, namely t — > ( of the form: 

d( = \(x,t)dt. (6) 

Such a transformation may distort directly the synchronization phenomenon acting on 
the times tt. To improve the timing between the oscillators given 7 ^> 1 and a < 1, we 
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perform the time transformation: 

, , s 7, if x\ 2 > and %2 i < and yi > 
I a, otherwise 

which shrinks the time between t\ and t\ enhancing a more accurate pacing between 
the oscillators. 7 may be chosen according to the pacing condition desired. For our 
purposes we fix 7 = 100. We can use the parameter a to control the average period. In 
the following we fix a = 0.11. The new time is given by (j = J* J A(x, t)dt. The equation 
of motion now takes the form 

£1,2 = A" 1 (x,t) [-0-1,2^1,2 - zi,2 + e(a?2,i - 2:1,2)], (8) 
£1,2 = A _1 (x,t)[ai i2 a;i,2 + 0.15yi j2 ], (9) 
i 1>2 =A- 1 (x,t)[0.2 + z lj2 (z 1 ^-10)] (10) 

The time transformation causes no changes in the state space, compare Figs. CD 
(a) and (b). However, the time series of x% x t and x± x ( are drastically modified 
[Fig. [T](b,e)]. Although, the time transformation is not able to interfere with the PS 
phenomenon [Fig. DDT)], it changes the frequency of the oscillators 



k)c = (<h)c « °- 998 



which implies that (Ti^)^ ~ 6.296. On the other hand, now max\({ — Q\ ~ (T^/420. 
Remembering that max\t\ — 1\\ ~ 2(Tj) t /3, we conclude that this time transformation 
yields an improvement of a factor of 280 for the timing. Of course, Eq. ([7j) can be altered 
to have an even better timing. These ideas can also be applied to a network. Whenever 
there is a cluster of oscillators in PS within the network, one can transform the time by 
Eq. (|SJ) suitably choosing A(x, t) to have a precise timing among all oscillators of the 
PS cluster. 



3. Phase Diffusion and Coherence 



Time transformation cannot destroy the synchronization. On the other hand, it does 
alter important characteristics of the dynamics, for example the coherence of the 
oscillators and the phase diffusion. By a time transformation we can transform an 
oscillator that originally is endowed with phase diffusion into an oscillator with an 
arbitrarily small phase diffusion. An interesting point is that in data analysis the phase 
diffusion plays a role in order to detect PS [19] . The idea is that one can detect PS by 
variations in the phase diffusion. 

In general the phase depends on the amplitude of the oscillator and the frequency 
can be written as: <fi(t) — u) + £(x, t), where uj is the average frequency of the oscillator 
and £(x, t), in many cases, acts as an effective noise due to the chaotic nature of the 
oscillator [8]. Therefore, the phase dynamics is generally diffusive, which means that for 
large time intervals one expects {[(pit) — wt} 2 )^ ~ Tt, where (-)^ denotes the ensemble 
average, and T the diffusion constant. Having the time of the ith cycle of the oscillator 
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Xj, we can write = i{Tj) + z/*, where (Tj) is the average period. By calculating the 
phase diffusion, we have 

- Uj(i(Tj) + v>)]\ = ^(H] 2 ) M « r,t 

Hence, t\ — i{Ti) gives the phase diffusion properties. 

Let us analyze the distortions in the phase diffusion by a time transformation and 
its effect on PS. Supposing that the oscillators xi and X2 are not in PS, we write 
t\ — t\ = a x % + £ l , where a,£ l G R are chosen to hold the equality. By performing 
a time coordinate change we endow the oscillator Xi with zero phase diffusion. This 
means that we have a new time ( with 

Aci = ci - cr 1 = i, (ii) 

i.e. AQ = At\/At\, where = t\ — t 1 ^ 1 . The new time coordinate is given by 



?1=0 

i 

71=0 

We have 

ici-GiHE( Af i- A ^)/ Af ii- ( 12 ) 

i 

Next, consider the maximum At\, namely 

maXiAt\ = 7 _1 . 

Thus, we have 

|£(Af'-Afi)/Afi|> 7 |f«-fi| (13) 

i 

which can be written as: 

ici-ai >7(«xi+r)- (14) 

Therefore, as the number of periods tends to infinity, the time event difference |Ci — C2I 
diverges. Thus, enhancing coherence in the oscillator does not introduce PS. 

4. Breaking down the Hypotheses on A 

By violating the conditions (ii) and (Hi), which guarantees the boundedness of A(x, t), 
PS can be introduced. Considering our former case where 

t[-i = axi + ?, (15) 

we could transform the time by 

A(x,t) = \ if t[ < t < t[ +1 . (16) 
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The timing condition is given by 



\C1-Q 




m— £ l 



A(x,t)| < \ai/i\ + \C/i 



(17) 



Thus, 



lim id-di<«, 



(18) 



■oo 



the time difference is bounded by a. Hence, this time transformation allows us to 
introduce PS between X! and x 2 . This does not contradict our results, because this 
time transformation is not bounded, violating the assumptions (ii) and (iii). When, 
one uses a non-bounded transformation A(x, t), like the latter one, the time is shrunk and 
becomes meaningless; there is no long term behavior with respect to (, since lim^oo Q 
is still bounded. The function A(x, t) can be made smooth without changing Eq. ([18]) . 
Here we have considered A(x, t) a step function without lost of generality. 

5. Feigning Phase Synchronization 

PS is invariant whenever A(x, t) fulfills the conditions (i — iii). However, it is important 
to mention that the phase definition must be defined consistently. This means that one 
must have the same phase definition before and after the time transformation, otherwise 
one could predict that PS is not invariant, due to the changing of the phase definition. 
Let us consider two spiking neurons Mi and A/a. We assume that the spike times t\ 
and f 2 are independent, with neuron Mi having a higher frequency, Fig. [2(a), in such 
a way that there is no n : m PS between Mi and A/2- Let t^ be the spike time of M2 
that precedes the i + 1th spike of Mi, Fig. Eta). Then, given Oi <C 7$, we perform the 
following transformation: 



This shrinks the time between t\ and t 7 ^ and stretches between t^ and creating 
bursts in A/2. In the rescaled time (, the bursts of M% is synchronized with the spikes of 
Mi. However, there is no synchronization between the spikes, since lim^oo |Ci — C2I ~^ 00 • 
On the other hand, after the time coordinate change, it is very tempting to introduce a 
phase for the bursts that increases 2n between two successive bursts of A/2. Changing the 
phase definition, the phase difference between Mi and A/2 becomes bounded. Therefore, 
it seems possible to introduce PS between two asynchronous neurons. However, this is 
a fake PS once that the phase definition is changed. 

6. Network Information Transmission 

Let us analyze the effect of time transformations in the information transmission 
in networks. For every pair of oscillators Xj and we can define a coordinate 




(19) 
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Figure 2. The effect of time transformation in two asynchronous neurons. In (a) 
the spikes in both neurons are governed by stochastic processes. They present no 
synchronization. From (a) to (b) we perform the transformation given by Eq. (119|) . 
which shrinks the time between t\ and t^* and stretches the time between t^* and . 
As a consequence, it seems that phase synchronization between the spikes in A/i and 
the bursts of Af 2 is enhanced. 



transformation 

x Sfc = x i + x fc ( 20 ) 
x i = x j - x fc, (21) 

that produces two positive conditional Lyapunov exponents (in units of bits/unit time) 
cr"(t) and cr- L (t). The mutual information rate (MIR) is bounded from above [12] 

Ic{t)<a\t)-a L {t) (22) 

The main goal is to know how the mutual information rate behaves as we implement 
a time transformation. By choosing a proper nonlinear A(x, t) in Eq. ([6]) we can 
introduce different time scales in the oscillators time series as well as endow the time 
transformation with as much information about the dynamics as we want. The main 
question is whether, under such nonlinear A(x, t), the information contained in A(x, t) 
could be transmitted to the oscillators. 

To answer this question we need to uncover the general transformation law for Iq- 
After some manipulations we can uncover the transformation law of Ic(t)' 

lc(0 < 7^> (23) 



where, again (-) t stands for the time average. For details see Appendix B 



Equation ( |23l) shows an invariant character of Ic- if Ic{t) = ; what happens in 
the absence of synchronization (correlation) between oscillators, no time transformation 
that respects conditions (i — in) can raise Ic(t) from zero. Hence, no matter how much 
information is contained in A(x, t), if there is no synchronization this information cannot 
be used. If, on the other hand, Ic(t) is positive, then Ic(() can be made arbitrarily 
large. 
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To illustrate our findings, we consider a network of four identical Hindmarsh-Rose 
chaotic neurons electrically coupled in an all-to-all topology, 

xj = yj + 3x 2 j - x| - zj + Ij + 22 Cjk{%k (24) 

k 

jjj = 1 - 5x| - y h (25) 

Zj = - rzj + 4r(xj + 1.6), (26) 

where Cjk stands for the coupling matrix. We use r=0.005, = 3.2, and random initial 
conditions. 




Figure 3. [Color online] Encoded time intervals between two spikes in neuron x 2 
for Cjk=0 (a), and for Cjk— 0.3 (b). In (c), the MIR between neurons Xi and x 2 , 
for the time-i (filled squares), and for the time-C (filles diamonds). BPS (burst phase 
synchronization) is found for Cjfc = [0.1,0.23]. In this regime only the burst are phase 
synchronized. PS is found for Cjfc = [0.23,0.25], and CS (complete synchronization) is 
found for Cjk — [0.25,0.3]. For Cjk £ [0.05,0.23] A(x, t) is smaller than 1, which provides 
an increasing in Ic up to 60%. For CS (A) f = 1.29 providing a decreasing in Ic- 



We define the following time transformation 

J a, if Xi = and y 1 > -4.6, 
A(Xl ' t) = \/?, if Xl = and yi < -4.6. (2?) 

This shrinks the time between spikes when the y\ > —4.6 and stretches when 
yi > —4.6 creating a frequency modulation between the spikes, which depends on the 
trajectory position. Hence, the transformation carries information about Xi. In our 
analysis we keep fix a = 0.5 and (3 = 2. t l j denotes the time of the ith crossing of the 
trajectory of Xj with the section Xj = (an spike event). The time interval between 
two crossings is At* = — t*. 

We introduce a symbolic dynamics which exhibits rather easily the results for the 
distinct synchronization regimes. We can encode the binary information about the 
transformation A(xi, t) by setting a to the symbol "0" and f3 to "1". Hence, we have for 
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two consecutive A = a: '00'; one A = a followed by A = (3: '01'; one A = (3 followed by 
A = a: '10'; and finally two consecutive X = (3: '11'. Whenever the time transformation 
is able to transmit the information about the symbols, we can access the information 
about x x in the spike time intervals of the other neurons. 

Figures [3](a-b) show return maps At l 2 vs. At 2 +1 of the neuron x 2 . We split this 
map into four return maps, depending on the value of the transformation A(x 1; t). That 
is, distinguished by the different symbols '00', '01', TO', '11'. The information about the 
values of A(xi,i) should be considered to be unknown, but here we make use of it to 
illustrate our ideas. 

By measuring At l 2 we should be able to infer the time interval At\, if the 
time transformation can transmit information. Figure [3](a) shows that return maps 
At 2 vs. At 2 +1 for the different values of A(xi,t) superimpose, and as a consequence it is 
impossible to discern whether the region that encodes for 00 is mapped to either 01 or 
00, and so on. That leads to a complete uncertainty about A(xi,£) by measuring At l 2 . 
Therefore, there is no exchange of information between xi and X2. 

The time scale of x 2 is being rescaled according to a function that contains 
information about the position of x 2 . From the way the function A is constructed, 
whenever y\ > —4.6 and X\ = 0, the oscillation frequency of the oscillator x 2 in the 
time-C frame is increased. Whenever y x < —4.6 and X\ = 0, the oscillation frequency 
of the oscillator x 2 in the time-£ frame is decreased. So, the oscillation frequency of 
x 2 is being modulated. Frequency modulation (FM) is a typical procedure to transmit 
information, a protocol in which the information signal is carried by the frequency of 
a wave. It would be natural to imagine that by modulating the oscillator x 2 using 
a time transformation based on the position of xi one could realize at least partially 
information about x x by making measures in x 2 . However, surprisingly, that is not 
the case in dynamical networks. Therefore, if elements in a dynamical network do 
not exchange information among themselves, there is no time transformation that can 
change this scenario. 

When the neurons are completely synchronized (for Cjk=0-3), we see in Fig. [3](b) 
that except for one point, the return maps At 2 vs. At 2 +1 for different values of A(x 1; t) 
are disjoint, which means that by measuring Af 2 we have complete knowledge about 
the trajectory of the neuron Xi. 

6.1. Effect of A( Xl ,t) onI c {t) 

We keep fix A(xi,t) and vary the coupling strength Cjk- Equation ( 1231) states that 
whenever (A(xi,t)) t < 1 the time transformation increases the MIR. In Fig. [3](c) we 
show the MIR between At\ and At 2 using the Shannon mutual information [20J, for the 
two time frames. Ic{t) denotes the MIR in the time-t frame and Ic(C) ^ ne MIR in the 
time-C frame. For CV,- e [0.5,0.23] (A(xi,t))t < 1 which provides an effective increasing 
in the MIR. 

In Eq. fT271) . A is defined to contain information about Xi. However, A could be 
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defined to contain information about an arbitrary information signal to be transmitted. 
In such a case, each disjoint region [as the ones shown in Fig. [3(b)] would encode 
information about this signal, which can be retrieved somewhere else in the network. 

A(x, t) can be constructed using information about some particular node of the 
network, a group of nodes. Whenever the oscillators are phase synchronized, we can 
improve the mutual information rate by using A(x, t) that contains information about 
the dynamics of the phase synchronized oscillators. 

7. Conclusions 

In summary, we have shown that for general dynamical oscillators it is neither possible 
to introduce nor to destroy PS by a time transformation. Furthermore, we have 
discussed possible application of these ideas to relevant technological problems such 
as nonlinear digital communication [IT] . Moreover, we have illustrated these results 
for nonsynchronized oscillators, showing that the enhancement of zero phase diffusion 
does not enhance PS. We have also discussed that breaking the boundedness condition 
imposed on A PS can be enhanced. However such a transformation is physically 
meaningless. Finally, we have shown that the time transformation can introduce the 
presence of distinct time scales, which can feign PS. Our findings might be relevant to 
several areas of natural science for the study of synchronization where the exact time 
the phenomenon took place is unknown and only a proxy for the time can be derived 
from the measurements. Examples can be found in geophysics when sediment cores are 
studied. Such situations may arise in dendrochronology, ice cores and three rings. 

Acknowledgment We would like to thank M. B. Reyes, M. Romano for a critical 
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Appendix A. PS invariance under time transformations 

We consider two general oscillators x,- = Fj(xj), where x 3 - G R™ and Fj : M n — > M. n , and 
analyze a general coupling scheme: 

xi, 2 = F 1)2 (x 1>2 ) + C 1)2 (x,t)[H 2il (x) - H 1)2 (x)], (A.l) 

where Hj(x) is the coupling vector function, and C^x, t) is the coupling matrix. Note 
that this scheme also takes unidirectional couplings (master-slave configuration) into 
account. We suppose that each x^ has a stable attractor and a frequency <pj = f2j(x, t), 
where f2j(x, t) is a continuous function (or Riemann integrable). Furthermore, we 
assume that there is a number M such that flj (x, t) < M. From now on, slightly 
abusing the notation we shall omit the dependence of the functions on the coordinates 
and on time, whenever there are no problems with the notation. Given a finite real 
number c, the condition for PS between xi and x 2 can then be written as: 

\Mt)-Mt)\<c (A.2) 
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First, we formalize the relation between PS and the timing condition given by Eq. (j5J). 
We have: 

I0i (t) - Mt)\ = I / 1 ftidt - [ 2 n 2 dt - / 1 n 2 dt + p\t)\ 

Jo Jo Jti 



where 



(3\t) 



Qidt 



tt 2 dt. 



Next, since </>j(tj) is equal to i x 2-n [2Tj . it yields: 

\<l>i(t) ~ Mt)\ < M\t\ - 4| + max0(t)\. 
The term maxj \(3 l \ is always bounded. By hypothesis we have 

.ill 



i*}-*h<a 



we have 



(A.3) 



< 



VLidt 



VLidt 





rt 




-.1 


n 2 dt 






i 








+ 


[ Q 2 dt 




Jt\ 



But now remembering that mctXjj {Qj{t)} = M, then 
1/3*1 < 2 MA 



(A.4) 
(A.5) 



(A.6) 



Therefore, a bounded time event difference \t\ —t\\ implies the boundedness of the 
phase difference. A similar argument shows that the boundedness of the phase difference 
implies the boundedness of the time event difference. Therefore, Eq. (jSJ) is equivalent 
to PS. 

We analyze the effect of time transformation PS. We assume A(x, t) to be (i) 
at least Riemann integrable (ii) finite, and (in) bounded away from zero. The two 
latter conditions are equivalent to the existence of two numbers 5 _1 ,?7 G M+ such that 
<5 _1 < A(x, t) < rj. Under the assumptions (i — Hi) we can demonstrate that PS is 
invariant under time transformations. First, we show that 

14-41 <«=H£-dl<« 

Noting that Q = J* 1 A(x, t)dt and Q = J <2 A(x, t)dt, we find 



ici-d 



A(x, t)dt 



A(x, t)dt 
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This may be written as \([ — Q\ — I Li X{x,t)dt\. However, since A(x, t) < r] we 

have 

ici-cii<^i4-4i- 

Thus, the boundedness of \t\ — t\\ implies the boundedness of \Q — Cll- Now, we show 
that 

ici-di<«=H*i-4i<«- 



We have 

\ L 1 L 2 



CI rQ 
A _1 (x,t)dt- / A _1 (x,t)dt 



which equals | j^ 1 A(x, t) l dt\. As 5 1 < A(x, t), we get 

if}-4i<*ici-&i 

Therefore, we conclude that there is PS in the "new" time frame ( if and only if 
there is PS in the original time t. 

The results stated in this section are general and do not depend on the attractor 
topology or coherent properties, as long as a phase can be introduced. Note that we do 
not have to know the phase equation, but only assume that it exists. 

The onset of phase synchronization, and even the phase equation, depends on the 
attractors topology and coherence [8]. If the attractor has a simple topology, that is, 
it has proper rotation, then, the onset of phase synchronization is given by a transition 
of the zero Lyapunov exponent to negative values. For such a case, the results of Ref. 
[T6] . concerning the invariance of the sign of the Lyapunov exponents, can be used to 
state the invariance of PS under time transformations. 

Appendix B. Transformation Law for Mutual Information Rate 

Representing one node dynamics of the network by 

§ = F(x), (B.1) 
the Lyapunov exponents of an invariant set of the phase space are defined as 

\yj\ 



hi = lim ln^, (B.2) 



t— >oo 

with 

yj = DFyi. (B.3) 

x(0) is a typical initial condition and y l tQ = yl(0)o~^ are the tangent vector at x(0). 

In Ref. [16j it is shown the transformation law for the Lyapunov spectrum. Under a 
time transformation A fulfilling the hypotheses (i — Hi) we have 

h < = A. (B.4, 
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Here, we show that the conditional Lyapunov exponents a» and a± follow the 
same transformation law, since the conditional Lyapunov exponents are the Lyapunov 
exponents of the network considering that all the initial conditions are equal, and 
therefore, the same results from [16] apply. 

Expanding Eq. (IB. II) linearly around the synchronous state s, and using the parallel 
coordinate defined in (120]) . we arrive that 

x 11 = 2[F(s) - DF(s)s] + DF(s)x" (B.5) 
= G(s,x"). (B.6) 

Proceeding in the same way for the perpendicular coordinate in f [2~Ti) . we arrive that 



x x = DF(s)x ± (B.7) 
= M(s,x ± ). (B.8) 

Now by means of G and M one can obtain the variational equations of Eqs. (IB. 61) 
and (IB. 81) . which provide the way small perturbations propagate along the parallel and 
perpendicular directions. From them, we obtain the Lyapunov conditional exponent 
erf along the parallel direction and the Lyapunov conditional exponent along the 
transversal direction, respectively. 

Then, by applying the results of Ref. [16], we conclude that 



(ai0) 

While the parallel conditional exponents are the Lyapunov exponents of the 
synchronization manifold, whose positive exponents measure the rate of information 
produced by the nodes if they were completely synchronous, the transversal conditional 
exponents are the Lyapunov exponents along the directions transversal (orthogonal) to 
the synchronization manifold, whose positive exponents measure the rate of information 
that can be erroneously transmitted between nodes. 

For a matter of simplicity in the notation, we denote the sum of all the positive 
parallel exponents by erf and the sum of all positive transversal exponents by . 

Then, by using the results of Ref. [12] , we can write 

Icif) < erf - ai, (B.ll) 

which in other words means that the mutual information rate, i.e. the rate with which 
information is exchanged between two nodes of the network, is given by the rate of 
information produced by the synchronous trajectories minus the rate of information 
produced by the desynchronous trajectories, the error in the transmission of information. 
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To intuitively understand Eq. (IB. lip , one can compare the right hand side of it with 
the usual definition of mutual information rate between a source of information, denoted 
by S, and the receiver of information, denoted by R, given by H(S) — H(S\R). The 
term H(S), which can be compared with erf, represents the rate with which information 
is produced in the source and the term H(S\R), which can be compared with , 
represents the rate of uncertainty remaining about the transmitted information after 
observing the received information, i.e., the rate with which information is erroneously 
transmitted. 

Then, taking into account Eqs. (IB. 91) and (IB.lOj) we have 

WQ < -jh-fc} ~ rt), (B.12) 
\ A )t 

concluding 

HO < j^- (B.13) 

[1] C. Schafer, M.G. Rosenblum, J. Kurths, and H.-H. Abel, Nature 392 239 (1998). 
[2] A. Riehle, S. Grim, M. Diesmann, A. Aertsen, Science 278, 1950 (1997) 
[3] R. Borisyuk, G Borisyuk, and Y. Kazanovich, Behavior and Brain Sc. 21, 833 (1998). 
[4] B. Blasius and L. Stone, Nature 406, 846 (2000); Nature 399, 354 (1999). 

[5] Steven H. Strogatz, Sync: The Emerging Science of Spontaneous Order, (Penguin Books Ltd, 
2004). 

[6] H. Fujisika and T. Yamada, Progr. Theoret. Phys. 69, 32 (1983). 
[7] L.M. Pccora, T.L. Carroll, GA. Johnson, ct al. Chaos 7, 520 (1997). 

[8] A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization: A Universal Concept in Nonlinear 
Sciences, (Cambridge University Press, 2001); S. Boccaletti, et al., Phys. Rep. 366, 1 (2002). 

[9] M. Le Van Quyen, et al., J. Neurosci. Met. Ill, 83 (2001). 
[10] M. Denker, et al, Neurocomp. 70, 2096 (2007). 
[11] T. Yalginkaya and Y.C. Lai, Phys. Rev. Lett. 79, 3885 (1997). 

[12] M. S. Baptista and J. Kurths, Phys. Rev. E 77, 026205 (2008); M. S. Baptista, J. X. de Carvalho, 
M. Hussein, "Optimal network topologies for information transmission in active networks", 
subm. for publication (|arXiv:0804.2983p . 

[13] P. Tass, M.G. Rosenblum, J.Weule, et al., Phys. Rev. Lett. 81, 3291 (1998). 

[14] F. Mormann, et al., Epilepsy Research 53, 173 (2003). 

[15] M. Marwan, M. Thiel, et al., Nonlin, Proc. Geophys. 9, 325 (2002). 

[16] A. Motter, Phys. Rev. Lett. 91, 231101 (2003). 

[17] B.R. Fayad, Ergod. Theor. Dyn. Syst. 22 437 (2002). 

[18] T. Pereira, M.S. Baptista and J. Kurths, Phys. Lett. A 362, 159 (2007). 

[19] Y.C. Lai, et al., Phys. Rev. E 73, 026214 (2006); M. Thiel, et al., Europhys. Lett. 75, 535 (2006). 

[20] We encode the spikes using the following rule. The z-th symbol of the encoding is by "1" if a spike is 
found in the time interval [iS, (i + l)8[, and "0" otherwise. A symbolic sequence is split into words 
of length L=4. The Shannon entropy is H = -^2 m P m log 2 P m , where P m is the probability of 
finding one of the 2 L words. The MIR is estimated by I c = [H(xi) + H(x 2 ) -H (xi;x 2 )]/(£ X L). 
We choose 8 £ [min(AiJ),max(Ai*)] to maximize Ic- 

[21] Note that this is also true for tj. This means that <j>j(t) = j Q j = 2n. Thus, in this construction, 
there is no contribution of the initial phase <pj (0) . 



